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In this work we study the local coupled Kuramoto model (LCKM) with periodic boundary con- 
ditions. Our main objective is to show how analytical solutions may be obtained from symmetry 
assumptions, and while we proceed on our endeavor we show some unexpected phenomena resulting 
from the symmetry properties: the existence of local attractors in the synchronized region; stabil- 
ity exchange for crossing fixed points; fixed point stability dependent on the manifold dimension; 
chaotic period and intermittent phase slips. As a result of our analysis, we show that stable fixed 
points in the synchronized region may be obtained with just a small amount of the existent solutions, 
and for a class of natural frequencies configuration we show analytical expressions for the critical 
synchronization coupling as a function of the number of oscillators, both exact and asymptotic. 
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I. INTRODUCTION 

Nonlinear systems tend to surprise us with their be- 
havior contrary to what we expect from knowledge, intu- 
ition and reasoning, even when the problem appears to 
be rather simple, as the Kuramoto model pQ and its local 
version. The applicability of the model for a variety of 
systems that range from physics [21 [3J , data mining [3] , 
neurosciences [5J |B] and others make it important to un- 
derstand the nature of unexpected results. Although the 
connection between the LCKM and applications is rea- 
sonably well established, the theoretical understanding 
is restricted to a few works [31 ITHT1] . Recently we stud- 
ied a locally coupled Kuramoto model (LCKM) above 
the synchronization transition [8]. We found a richness 
in the solutions which was not expected, mostly because 
intuition tells us that a ring of oscillators should behave 
as a chain for large systems. We studied how they are 
born and the region of phase space where the solutions 
exist. The complexity of the problem makes it extremely 
hard to analyze therefore we shall start looking at small 
systems with increasing difficulty to see if we can obtain 
some results from which we can infer large N behavior. 

Our objective is to calculate the behavior of the solu- 
tions for several cases where the distribution of frequen- 
cies have well determined symmetries. We will study how 
these solutions appear and follow them as the system size 
increases. To do that we shall start with the smallest sys- 
tem (N=3) which has already received the attention of 
Maistrenko et al [151 ES] an d we are going to increase 
the size in the search of general properties. Obviously we 
shall not attempt to have an encyclopedic coverage but 
will show the difficulties encountered in solving the gen- 
eral problem with any arbitrary distribution. In doing 
so we have come across some results on phase slip which 
happens at the transition to full synchronization which 
seem to contradict previous results for specific symme- 



tries [TUl [HI , which we shall discuss along the text. 

The manuscript is organized as follows: the model is 
introduced in Section II; in Section III A we present a well 
known case, that of three oscillators, where we also obtain 
for the first time the analytic result for the dependence 
of the solutions on the coupling constant above synchro- 
nization; after that, in Section III B, we consider a sym- 
metric N — 4 case and show a general method for solving 
the equations for any number of oscillators, along with 
its difficulties; in section IV we implement the method in 
a very organized system, composed of only two natural 
frequencies, and show some unexpected features concern- 
ing the dynamics and stability; the analyticity limits of 
the model are explored in section V within the context 
of a specular symmetry, and we show how the model's 
degeneracy may be lifted when the symmetry is broken; 
in section VI we show some examples on how organiza- 
tion on the natural frequencies may lead to analytical 
(asymptotic) expressions for the critical synchronization 
coupling; a summary of the results and some possible 
further extensions are left for the final section. 

Before ending the introduction we would like to define 
some of the notation used throughout the text. When 
treating the bifurcations of the model we are not inter- 
ested on their classification, but on the characterization 
of the fixed point's stability. With that in mind, we used 
a simplified notation: the bifurcations are named accord- 
ing to the type of fixed points created at them; for the 
fixed point classification, we called a saddle every solu- 
tion containing both positive and negative eigenvalues 
of the Jacobian matrix and stable (unstable) node when 
every eigenvalue is negative (positive), regardless of the 
imaginary parts. Another important definition concerns 
the inverse sine function: unless specified on the text, 
arcsin (x) is to be considered as defined solely on the in- 
terval [— 7T, 71"]. 
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II. PERIODIC BOUNDARY CONDITIONS: 
EFFECTS OF SYMMETRY 



some of the problems within this description when ex- 
tended to higher number of oscillators. 



The LCKM under periodic boundary conditions 
presents a rich landscape of solutions as discussed in ref- 
erence [S] . From now on we shall follow the nomenclature 
used in that reference. The system is described by the 
set of equations: 



AT [sin ( 



?n+l 



,)], (1) 



for n = 1, ...,N, where u) n 6 {^} ^ IS the set of natural 
frequencies. The topology of the ring is defined by the 
conditions 8m+i = Q\ and da = 6m- Similarly to the case 
of the chain |17j there is a minimum value of the cou- 
pling constant K for which the system synchronizes to a 
common frequency f2 = jv T^-i Under full synchro- 
nization the set of equations p| can be written as 



K 



sin i 



smc 



1. 



,N, 



(2) 



where 0„ — 9 n — 9 n+ \. The condition for the phase 
differences 0„ (AT, {oj}) to lock depends on the number 
of oscillators N and the value of the coupling constant 
in the region AT > K s . In the synchronized state every 
variable <p n can be written in terms of one n « arbitrarily 
chosen: 



sin < 



1 " 



for n = l, 



— 1 and 



sin m n = sm i 



1 ™ 



(3a) 



(3b) 



j=n*+l 



for n = n* + 1, AT — 1. Since the identity X^=i <Aj = 
allows us to write 0jv as a sum of all the other phases n , 
with n = 1, ...,N — 1, the set of equations ([2]) is reduced 
to a single equation on two variables ((f) n *,K): 



sin 




sm < 



A" 



(4) 



The description of the solutions as a function of K is not 
a simple extension of that of the chain since for the ring 
topology the equations for the phase differences are not 
independent. It is necessary to analyze simple particular 
cases where the small size, or particular symmetries, will 
give us a hint on how to proceed to the general case of 
randomly selected natural frequencies. This is what we 
shall do in the following sections. 



III. SMALL SYSTEMS ANALYSIS 

In this section we make an analysis of small systems, 
particularly the N = 3 and N = 4. We show two ap- 
proaches for obtaining analytical solutions, and discuss 



A. First non trivial case: N=3 

The first non trivial case is a system of three oscillators 
and natural frequencies oj\ — — L03 — uj and u>2 = 0. 
Clearly O = and the synchronized region is described 
by: 



4> 2 ) + sin 0i 
n 4>2 — sin <f>i 



0. 



Since (5b) has only two solutions, 



(5a) 
(5b) 

(6) 



it is possible to use only one phase difference 4>\ = </> and 
rewrite ^ as a single equation K ((/)) for each solution: 

Rl (*) = • ^ > K " w - —a, ■ ( ? ) 

sm (p + sm 20 sm (p 

The function K 11 ((f)) has the same structure as a chain 
of oscillators, with just a single minimum at K 11 = oj 
with sin0 = 1. On the other hand K 1 (0) has two min- 
ima for K > 0. If we define z — sin</> these minima can 
be written as K± — f (z±), i.e.: 



4 (ii 2 ^ 



(8) 



where 



15 



'15 



33 



32 



32 



(9) 



Each of the three minima gives birth to a pair of solu- 
tions for K > K s , and their stability is characterized by 
the eigenvalues of the Jacobian matrix. When we write 
the equations of motion as a function of the phase differ- 
ences 



- K [2 sin 0! 



— sin 



yi — oj — K [2 sin 0i — sin 02 + sin ((fix + <j) 2 )\ , (10a) 
02 = + AT [sin 0i — 2 sin 02 — sin (0i + 02)] , (10b) 

the eigenvalues for each type of solution are defined by 

A±(0) = K [-2 cos - cos 20 ± (cos 0- cos 20)] ,(11) 

\± (0) = — (2 ± V2 y/3 - 2 cos + 3 cos 20) . (12) 

To analyze the stability of the system, we plot the func- 
tions 1/K±(4>) and 1/K 11 (0), with the corresponding 
eigenvalues of the Jacobian matrix versus in figure [T] 
We notice that the first maximum at represents a 
saddle- (stable) node bifurcation: two solutions are born, 
one stable with 0—^0 and one unstable with —> 2tt/3 
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fifi K i = Ks (a) /"(<*) K n (b) 




FIG. 1: (Color online) (a) Plot of the functions 1/Ki (0) 
where the stable solutions are shown in solid black lines and 
the unstable solutions in dashed lines and \± (0) (thin blue 
and red lines) for ui — 1. Saddle-(stable)node bifurcation at 
K\ = K a , the solutions are stable on the left of the bifur- 
cation and one of the solutions is unstable, on the right. At 
the minimum of KL there is a saddle-(unstable)node bifurca- 
tion, (b) For 1/K 11 {(f)) and A± (0) we see a (unstable)node - 
(unstable)node bifurcation at the minimum of K" (0), which 
become saddles at the edges. 




FIG. 2: (Color online) General solution K± (02) represent- 
ing the fixed points of the N — 3 system. Stable (unstable) 
branches are represented by solid (dashed) lines, (a) Com- 
plete solution for the symmetric case shown in figure [l] where 
it is possible to see a pitchfork bifurcation at 02 = 7r/2. (b) 
General structure of the synchronized region for a configura- 
tion of natural frequencies without symmetries: wi = 1 and 
UJ2 = 1/2. The six fixed points solutions come from three 
different bifurcations. 



in the limit K — > oo. The solutions generated at the 
minima KL (figure |TJl) and K 11 (figure [T]d) are unsta- 
ble, since the real part of \{_ and \+ is always positive: 
K l _ is a saddle-(unstable)node bifurcation and K 11 is a 
(unstable)node-(unstable)node bifurcation. 

For this system with just a few oscillators it is not 
necessary to impose symmetry properties on the natural 
frequencies in order to obtain an analytical description 
of the synchronized region. If we consider the natural 
frequencies to be randomly chosen, it is always possible to 
assign positive values to uj\ and ui 2 and set W3 = — u)\ — u> 2 
(for S1 = 0). With this configuration the fixed point 
solutions are described by the equations: 



sin (0i + 2 ) + sin (/) 2 = 
sin 0i — sin 02 = 



. + UJ 2 

K : 

L0 2 

Jc' 



(13a) 
(13b) 



A suitable manipulation of equations ( 13 1 enable us to 
write a single equation for 02, 



sin 



1 - (sin 02 - uj 2 x) 



[(ui + uj 2 ) x - sin 02 - (sin 02 - uj 2 x) cos0 2 ] 2 , (14) 



where x = K . Since (14 1 is a polynomial of second 



order in x, we find solutions in the same fashion as equa- 
tion Q, i.e., functional forms K~ x (0 2 ) describing all the 
fixed points on the synchronized region: 



Kg 1 



> 2 ) = 

(wi + 2w 2 ) sin 02 (1 + cos0 2 ) ± \A (02) 

wf + 2co 2 (wi + lo 2 ) (1 + cos 2 ) 
>a) = sin 2 02 [w 2 cos 2 2 

+2cj 2 {0J1 +uj 2 )(1 + cos 2 )] . 



(15a) 



(15b) 



For the symmetric case the fixed points are those shown 
in figure [T] 



It is now possible to see that the K 11 solutions are born 
from the rightmost branch of K 1 (figure [2^l) , which lead 
to the pitchfork bifurcation description by Maistrenko et 
al. |15) . A general configuration of natural frequencies 

f aerates the structure of fixed points shown in figure 
, where the maxima of -fCjl 1 (02) appear separate, thus 
aracterizing the presence of three bifurcations on the 
phase space. Unfortunately, the analytic expression for 
K s {uj\,u)2) does not admit a simple representation and 
to best of our knowledge it has not yet been reported. 



B. N=4 with specular symmetry 

Based on the same assumptions as the previous case, 
we shall start analyzing the ring with N = 4 oscillators 
where the configuration of frequencies follow a prescribed 
symmetry. It is woth mentioning that Maistrenko et al. 
[Td] have already studied the case of N=4 for the full 
connected Kuramoto model. We treat here the case when 
the oscillators present specular symmetry on the natural 
frequencies, with W3 = — lo 2 , — and both oj\ and 
o>2 positives. 

The symmetry in the equations imposes sin0i = 
sin 03 , with solutions 



(16) 

For the 0^ case, the equations in the synchronized state, 



sin (n + 02 ) + sin 0i = 



K 1 



sin 02 — sm 0i 



oj 2 
~K 



, (17) 



permit the existence of a solution only if lj 2 = — u)\ . As 
this case belongs to a more general scenario described on 
the next section, we will consider only the 0| solution, 
whose synchronized region is described by the equations: 



sin (20i + 02 ) + sin 02 = ^ t— ^ , (18a) 



K 



sm 02 — sm 0i 



uj 2 



(18b) 
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After some manipulation, equation (18a I can be written 
as 



sin 



!>i (1 - sin^ 



UJ-2 



2K 



(1 - sin z : 
- sin 02 (l 



sm 



(19) 



If we would like to search for a solution of the form 
(</>2), as it was made in the case N = 3, equation 
( 19 1 shows that it is necessary to obtain the roots of a 



vestigation shows that if we define z = sin 02, it is possi- 
ble to write a third order polynomial, 

^ (y2 - Vi) - z 2 y 2 (3y 2 - 2y x ) + z (3y| - y x y\ +y%- y 2 ) 
1 



(% 2 4 



3y 2 + y{ + 2yiy 2 = 0, 



(20) 



where y x . 2 — ^12/ K. On this slightly different approach, 
the fixed point solutions come as the roots of equation 



fourth order polynomial. However, a more accurate in- (20), namely, 



ri{K) = 



2<jJ\U) 2 



3w| 4w 2 wf + 12K 2 (uii 



+ if 4 /o /3 (K) 



3K (u>i — uj 2 ) 



6X3 ( Wl _ W2 )/ 1 /3 {R) 



2uj^ 2 -3u? 2 -4(l±iV3) 3X 2 K-^) 2 +^ 2 2 +{-l±iV3K*)f 2/3 (K) 



3K (wi — uj 2 ) 
where the function f (K) is defined by 



12^3 ( Wl -W2)/ 1/3 (K) 



fo(K) = fx (K) + 3\/3V/ 2 (K), 



h{K) 
h{K) 



ujiUJ 2 (<jJ\ — oj 2 ) 2 (3oj 2 — 2bj\uj 2 + 3w|) 



'6 



9 



16 



K 

Ulfu! 2 (0J\ — UJ 2 ) 4 



K 4 

(cji - lu 2 ) 4 (27wf - 36w 3 w 2 + 2u]\u]\ - 36wio;| + 27w|) 



r 



(21a) 
(21b) 

(22a) 
(22b) 



m^L. (22c) 



The roots r x (K) and (K) give the locked solutions 
of sin 02 (K) for K > K s , but as | r\ | is always greater 
than 1 this root is not a valid representation of sin 02. 
From these solutions it is possible to build a bifurcation 
diagram for the system once we know the stability of the 
solutions. Bearing in mind that we are searching for sim- 
ple analytical expressions, we will consider a reduction of 
the phase space dynamics. With a proper set of initial 
conditions, i.e., 3 (t = 0) = 0i (t = 0), with X (t = 0) 
and 02 (t = 0) randomly generated, the dynamics is set 
to occur in a 2-dimensional symmetric manifold. Now, if 
we follow the same procedure of the previous section and 
write the equations of motion for the phase differences 



01 = OJl — UJ 2 

-K [2 sin 0i - sin 2 + sin (20i + <j> 2 )} , (23a) 

02 = 2lo 2 + 2 (sin 0i - sin0 2 ) , (23b) 

we obtain a closed analytical expression for the eigenval- 
ues of the Jacobian matrix. 



A+ = 



t±Vt 2 - 4A 



(24) 



where the trace r and the determinant A of the matrix 



are given by 

r = -2ifcos0i [1 + 2 cos (0i + 2 )] , (25a) 



A = -4X 2 cos 2 f ) [1- cos 20i 



-2 COS (01 +02)] 



(25b) 



At this point we find a problem: the locked solutions 
are given in terms of sin 2 while the eigenvalues of the 
Jacobian depend on 0i and 02- A way out of this conflict 



can be reached if one notices that for each root in (21) 
there are two solutions for 2 : 



02 (ri) = arcsinri, 
02~ (r 2 ) = arcsin 



b T 2 (ri) = 7r — arcsinri, 



(r 2 ) = 7T 



arcsm r 



(26a) 
(26b) 



do not depend on 0i, it is necessary to go 



Since r\ < 

back to equations ( 18 ) to determine the dependency of 
0i on 02 for each solution of ( 26 1 . A numerical analy- 
sis of the equations shows that 0i = arcsin (sin 02 — ^) 
leads to the solutions 02 (ri) and 02 (r 2 ), while 0i = 
7r — arcsin (sin 02 — ^) corresponds to the complemen- 
tary branches 0| J (ri) and 0| J (r 2 )- 

As a result of this analysis it is possible to observe 
that, while the graph sin 02 x K shows K s as a bifurca- 
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FIG. 3: (Color online) Roots n (K) and (K) of (20 1 rep 



resenting the bifurcation diagram in the synchronized region 
for lui = 1 and tJ2 = 1/3. Full (dashed) lines represent stable 
(unstable) solutions, a) Bifurcation diagram projected on the 
space sin 02 x K showing the bifurcation at K a where two 
solutions are born: one stable with sin 02 — > and another 
unstable with sin 02 — > 1 in the limit K — > oo. b) Bifurcation 
diagram showing the simultaneous birth of the solutions in 
different regions of phase space at K s , a saddle-(stable)node 
in the region cos 02 > and an saddle-(unstable)node for 
COS 02 < 0. 




FIG. 4: (Color online) Phase space showing the evolution of 
the fixed points as a function of K for a fixed configuration 
of natural frequencies. Black(thick) and blue(thin) lines cor- 
respond to 0i = (02 = 0). a) K = 0.7: just before K s , 
we can see the simultaneous bifurcations, b) K = 1.5: stable 
node and saddle created in the region with cos 02 > (cen- 
tral part of the figure), and saddle and unstable node created 
in the region cos 02 < 0, in the upper part of the figure, c) 
K = 20.0: for large values of K the stable node moves to- 
wards (0i,02) = (0,0), the two unstable saddles go towards 
(7r/2,7r/2) to merge into a unique saddle, d) In the limit 
K — > oo, the unstable node moves towards (ir, n) while the 
curves deform to be a saddle at (0, 7r), when K — > oo. 



tion giving birth to two solutions (figure [3^,), the identity 
sin0 = sin (tt — 0) implies that for each branch of sin 02 
there are actually two fixed points in the 2 x K space 
(figure [3J)). The unexpected (though straightforward) 
conclusion is the feature that the K s bifurcation occurs 
simultaneously in two regions of the phase space, but the 
fixed points do not share the stability properties: in the 
region defined by cos 02 > it bifurcates into one sta- 
ble 2 (r%) and one unstable branch 02 (f 2 ), meanwhile 
for cos 02 < it has two unstable solutions, a saddle 
(r"2 ) and an unstable node [r-y). 



Now it becomes evident why we chose a reduced sub- 
space for the dynamics: due to its lower dimensional- 
ity it is possible to visualize the evolution of the null- 
clines and the creation of the fixed points in the phase 
space, as shown in figure [4] Just before the bifurca- 
tion at K s it is possible to visualize the two simultane- 
ous fixed points being formed at (0i,02) ~ (0.4, 1) and 
(0!,0 2 ) « (2.5,2.5) (figure^). Figured (K > K.) 
shows the stable node at (0i, 02) ~ (0, 0.5) and the sad- 
dle at (0i, 02) ~ (1, 1.4) born in the region cos 02 > 0, as 
well as the saddle at (0i, 02) ~ (2.5, 1.5) and the unsta- 
ble node at (0i, 02) ~ (3, 2.5). In the limit K — > 00 the 
stable node shifts towards (0,0), the two saddles move 
towards (tt/2,%/2) until they merge into a unique sad- 
dle, and the unstable node shifts in the direction of (it, it) 
(figures^ and|IJi). 

For this specific system the reduced manifold contains 
all the solutions of the complete 3-dimensional phase 
space, and the stability of the fixed points does not suffer 
any alteration when one reduces the dimension of the dy- 
namics. However this is not the general behavior for the 
LCKM when symmetries on the natural frequencies are 
taking into account, as we will show on the next section. 

The value of the coupling constant at the synchroniza- 
tion transition can be obtained as the minimum of the 
function K (sin 02) = K (z), such that d z K (z) \ z=z * = 
implicitly calculated in ( 20 ) gives two solutions 



uj 2 (2u>i — 3^2) ± \J + 3K 2 (uji — Lo 2 y 



3K (wi — co 2 ) 



(27) 



Once we take into account the two roots z* ± in (20 1 we 
obtain four solutions: 



K±> ± (u Jl ,u J2 ) = ±—^ Jk±("i,"2), (28) 

»V2 \u)i - uj 2 \ v 



where the function (uii,(jj 2 ) is defined as 
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(u)i,u>2) = 27uj\ — %uj\ (au>2 ± \J k\ (wi, w 2 )^ + U2Lo\ {2u>2 ± 5y/k~i (wi,u; 2 )^ — ^36w 2 =F (wi, cl> 2 )^ 

+9wf (3oj 2 T Vki (wi,w 2 )) , (29a) 
ki (wi, w 2 ) = 9cJi - 14wiw 2 + 9cj 2 . (29b) 



For the case considered (uji and w 2 positive) the solutions 
Kf ,+ are always complex conjugate pure imaginary num- 
bers, and do not represent a physical solution. On the 
other hand the solutions Kf'~ are always real, but, since 
K~'~ is always negative the critical coupling at synchro- 
nization is given by 



k« = k: 



(wi,w 2 ) 



(30) 



In the figures we have taken u)\ = 1 and cj 2 = 1/3 for 
which K s « 0.734. 

With this two examples of small systems (N = 3 and 
N = 4) we gave a description of how analytical expres- 
sions for the location of the fixed points and the criti- 
cal synchronization coupling may be obtained, which we 
summarize here: starting from equation Q, if we make 
the trigonometric expansion of the sin^jv term into sines 
and cosines of all the phase differences, then by squaring 
the two sides of the equation (with subtle rearrangements 
of the terms) it is possible to use equations ^ to write a 
polynomial equation either for z = sin0„* or x = K . 
Ultimately the fixed points will come as the roots of the 
polynomial, and the critical synchronization coupling K s 
appears as the first minimum of the function K (</>n«), 
computed with respect to either <j) n * or sin </>„* , as long as 
</>„» is chosen properly (for further information about this 
point see [5]). The problem with this approach is that 
the degree of the polynomial increases with the number 
of oscillators, and for N > 5 the degree is typically too 
high to obtain analytical roots. Nevertheless if we make 
use of highly symmetrical configurations of natural fre- 
quencies the behavior of the system may be analytically 
explored to a large extent. This is the task that we shall 
undertake in the next section. 



IV. TWO NATURAL FREQUENCIES 

The simplest system with nontrivial behavior one may 
consider is composed of only two natural frequencies. But 
as local coupled systems present a strong dependence on 
the local configurations, we will get rid of the inhomo- 
geneities by assuming an organized distribution of the 
natural frequencies: 



UJ. n 



uj, for n odd 
— uj, for n even 



(31) 



In the synchronized region the symmetry reduces the 
equations to a set of constraints on the phase differences 



sin < 



sin (pi, 
sin 4> 2 , 



for n 
for n 



3,5,. 
4,6,. 



while 4>\ and </> 2 satisfy 



' N-l 

E < 

.n— odd 



N-2 

E 



sin (pi 



K 



(32) 



,(33a) 



sin <j>x — sin <f>2 — —.(33b) 
K 



A straightforward solution of this set corresponds to take 
4> n = 4>\ for n even, <p n = </> 2 for n odd and (f>% — —<fii, 



such that ( 33 ) can be written as 



sin < 



2K' 



(34) 



As this set of solutions (4> n ) maximizes the left hand side 
of (33a), the critical synchronization coupling is deter- 



mined as the first solution of equation (34), i.e 



2 



(35) 



Although this specific choice of the phases gives a fixed 
point solution valid for K > K s , the synchronized re- 
gion is full with solutions obtained from combinations of 
4>ti = 0i,2 and 4>~ — 7r — 0i,2 ; in a way that a complete 
description of the system involves all possible combina- 
tions of the solutions to equation (32). 



To build the space of solutions, we will explore the 
system's symmetries to get the simplest description of 
the dynamics. Since equation ( [32] ) tells us that there 
are two independent phases, it is possible to consider 
proper initial conditions <p n (t = 0) that reduce the origi- 
nal (N — l)-dimensional system down to a 2-dimensional 
set of equations. If we define 



cf>i, for n = 3, 5, N — 1, 
2) for 71 = 4,6,..., iV- 2. 

7r — <j)\, for 71 = 3, 5, N- 
7r — (p2, for 71 = 4, 6, N- 



(36a) 
(36b) 



the equations of motion are written as 

4>i = 2w — K (2 sin <f>\ — sin <p2 + sin , (37a) 



-2uj + 2A" (sin (j>\ — sin 2 ) 



(37b) 
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where 



= 91+92 



N-l 

£ 

n=3,5,. 
■)N~3 



N-2 



^(0 2 ), (38) 



7i=4,6,... 



defines each of the 2 6 solutions that compose the 0at 
term, thus generating all the subspaces of the system. 
The stability of the solutions is described by the eigen- 
values of the Jacobian matrix 



A, 



4A 



(39) 



where r, the trace and A, the determinant are defined 
by 

t = -K [2 (cos 4>i + cos 2 ) + cos <i>(90 1 $] , (40a) 
A = 2K 2 [cos <& (cos 0290!$ - cos fad^®) 

+ cos 0i cos 2 ] • (40b) 

Starting with a N = 6 system, the symmetrical sub- 
spaces are defined by <£> ±±=t such that the equation (33a 
has the following representations in each subspace 



$+++ 

$-++ 
$++- 

$+— 

$-+- 

$— + 
$+— 
$ 



sin (30i + 202) + sin 0i 



K 



(41a) 



- sin (0i + 20 2 ) + sin 0i = — , (41b) 

K 



— sin 30i + sin 0i = — , 
K 

sin ( -0i + 20 2 ) + sin 0i 



K ' 



sm 0i 



2K' 



(41c) 
(41d) 

(41e) 



Since the set (41) corresponds to all possible solutions 
of equation (33a I, every fixed point in the synchronized 
region can be represented by the functions K^ 1 (0i), ob- 
tained from pi]) and (33b). The solutions, with their 



corresponding bifurcations point in each subspace repre- 
sented by the pair (K* , 0*) , are described in the table [I] 
and illustrated in figure [5] with u> = 1. 

An unexpected feature that comes from the symmetry 
presented by the system concerns the characterization 
of the critical coupling: our analysis shows that K s is 
highly degenerate, since each subspace has two bifurca- 
tions responsible for the synchronization: a node-node 
bifurcation (with opposite stabilities) coming from the 
< + solutions, always accompanied by a saddle-saddle or 
a saddle- (unstable) node coming from the 0^ solutions. 
It is important to notice that not only at K s but at all 
values of the coupling constant where new solutions are 
born there exists the same composition of double bifur- 
cations, as can be observed in figure [5] 

Although the solutions found on the subspaces corre- 
spond to all the fixed points of the system, the stabil- 
ity obtained from the eigenvalues of the bidimensional 



Subspace $ +++ 


vi = if ni - fa 


711 = 
71l — 1 
711 = —1 


-> (if ,0*) = (w/2,ir/2) 
-+ (if ,01) = (w,57r/6) 
-> (if ,0*) = (w,ir/6) 


02 = 7T (2ri2 — 1) 
-301 


712 = -> < 


(if ,0i) = K2 )7 r/2) 
(if,0!)« (1.84,3.56) 
k (if ,0*)« (1.84, -0.42) 


Subspaces $ ++ and $ ++ 


<f>2 = 27T71i — 01 


Til = 0- 


+ (if ,0*) = ( W /2,7r/2) 


A - _ 7r(2n 2 + l)-01 
92 — 3 


712 = - 
712 = 1 " 
712 = 2 - 


^(if,0*)« (1.84, 1.88) 
^(if,0t)« (1.84, 1.26) 
■+(jr,«) = (w/2 ) 7r/2) 


Subspace $ + + 


02~ = 301 

02~ = 7T — 301 


(if 

(if 
(if 


,#) = (w/2,jt/2) 
,0*)« (1.84, -0.42) 
,0*)« (1.84, 3.56) 


Subspace $ + 


j+ _ 27rn 1 +0i 
92 — 3 


n 1 =0 

71l = 1 
Til = — 1 


-> (if*, 0i)« (1.84, 1.26) 
-»• (if*, 01) « (1.84,1.88) 
->(/r,«) = (w/2,ir/2) 


02 = T (2712 ~ 1) 
+01 


n 2 = 0- 


-> (If ,01) = (w/2,ir/2) 


Subspaces $ + , $ 


H and 


<t>2 = -fa 

0J = -7T + 01 


= (w/2,7r/2) 



TABLE I: Solutions for each symmetric subspace and the re- 
spective bifurcation points. The solutions 0^ (0i) for each 
subspace (on the left) depend on integer numbers 711,2- When 
the two solutions have the same bifurcation points (if ,0*), 
the subspace is completely symmetric with respect to 0^ (ex- 
cept for the stability), and the number of bifurcations is dou- 
bled. 



Jacobian matrix in equation ( 39 ) does not necessarily 
correspond to the stability of the solutions in the com- 
plete phase space. Thus, in order to know the complete 
stability of solutions it will be necessary to analyze the 
eigenvalues from the (N — 1) x (N — 1) Jacobian obtained 
from all the equations of motion instead of the bidimen- 
sional set. 

From the analysis of the fixed points born at K s we 
observe the existence of only two stable solutions, both 
of them originated from $ +++ : the stable solution from 
+ (ni = 0) and one of the saddles from 9"^ ( m this case 
the perturbations to calculate the eigenvalues are re- 
stricted to the symmetric subspace). The latter is born 
at a saddle-saddle bifurcation in the subspace (figure [6^,), 
but it gains stability in the complete phase space as the 
bifurcation becomes a node-node type (figure ^jp shows 
the eigenvalues obtained numerically for the two solu- 
tions, with K slightly above K s ). All other fixed points 
born at K s were created from bifurcations saddle-saddle 
like in the complete phase space. 
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FIG. 7: Time evolution of the instantaneous frequencies 9\ (t) 
and 82 (t) just below synchronization (K = 0.499). (a)For a 
set of initial conditions seemingly chaotic behavior of the time 
difference between slips is observed, (b) Time evolution of the 
instantaneous frequencies for a set of oscillators which wander 
around phase space and only eventually find a neighborhood 
of a bifurcation. 



FIG. 5: (Color online) Bifurcation diagram K~ [<f>{) show- 
ing the solutions and their stability defined in the symmetric 
spaces. Solid (dashed) lines represent stable (unstable) solu- 
tions, (a) Subspace <I >+++ . (b) Subspace <E>~ ++ . (c) Solutions 
<f>2 of subspace "J?" 1 . (d) Subspace $ 1 . 
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FIG. 6: (Color online) (a): Solutions fa (<j> x ) of $+++ born 
at K s through a saddle-saddle bifurcation in a bidimensional 
subspace (dashed lines). Eigenvalues of the Jacobian matrix 
(thin lines), (b) Complete set of eigenvalues of the Jacobian 
matrix for the whole phase space, at the bifurcation the sta- 
bility of the saddle changes into a node-node, the solution for 
4>i — > 7r/4 (as K — ¥ oo) is a stable node (black dots), while 
the solution cj>i — > 3n/4 is an unstable node (blue squares). 



The fact that a stable solution in a given subspace loses 
its stability in the complete space is not a surprise, since 
the subspaces may represent only the stable manifolds of 
a saddle. On the other hand the fact that a fixed point is 
unstable in a given subspace and stable in the complete 
space is not yet understood. 

The stability analysis in the region slightly above K a 
gives us unexpected results: from a total of 32 solu- 
tions created from 16 different bifurcations (remember 
that there exist 8 subspaces of symmetric solutions, each 
with two bifurcations), there are two stable fixed points 
originated from node-node bifurcations while the other 
28 are saddles. If we consider a system with such a large 
number of saddles and two stable nodes, it will be only 
natural to expect an unusual behavior. We analyzed the 
time evolution of the instantaneous frequencies which is 
illustrated in figure [7] When synchronization is about to 
happen (K = 0.499) we expected to observe phase slips, 
with the period between slips proportional to (K a — K) 



[TTJl [T7] . Surprisingly what we observed here is a strong 
dependence on the initial conditions: from some initial 
conditions the trajectories of the oscillators pass near the 
regions where the bifurcations appear, but as there are 
several possible bifurcations they get shifting from one to 
another (each producing a slip of a different strength) and 
the period between phase slips behaves in an apparently 
chaotic way (figure [7^l) ; at the same time there are ini- 
tial conditions that make the oscillators wander around 
phase space and only eventually find a neighborhood of 
a bifurcation, what does not necessarily happen simulta- 
neously for all oscillators, resulting on erratic behavior 
and intermittent phase slips as seen in figure [TJa. 

This unusual behavior is also observed for values of K 
slightly above K s : with a given a set of initial conditions 
9 n (t = 0) = 77, where r\ is a random variable denned from 
a uniform distribution in the interval [—a, a], for values 
of a less than 1 (approximately) the system always finds 
the stable fixed point of fa (n\ = 0); for a > 1 there ex- 
ist initial conditions that do not drive the system into a 
synchronized state, where all the oscillators remain con- 
stantly searching for the fixed points (figure [tJd shows the 
typical behavior of the instantaneous frequencies in this 
case); the fixed point of fa has an extremely small basin 
of attraction, very near initial conditions are necessary 
for the fixed point to be found. This behavior implies 
that the stable nodes born at K s correspond only to lo- 
cal attractors in the phase space. A possible explanation 
for this phenomena is that the surfaces defined by fa = 
form a cocoon around these fixed points, acting as a re- 
pulsive region for the trajectories started outside. On the 
other hand the high dimensionality of the system makes 
the task of verifying this hypothesis herculean and out- 
side the scope of this manuscript. 

So far we have described the unusual behavior of the 
system prior to and on a neighborhood immediately af- 
ter the onset of synchronization. However a set of unex- 
pected features also occurs on the solutions in the syn- 
chronized region, whose peculiarities come from the high 
density of solutions. In order to show some of these fea- 
tures we shall focus our attention on the $++■■• subspace, 
which present most (if not all) of the features and has a 
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direct link to the system's stability in the whole phase 
space which will be explained in the next section. 

For an arbitrary number N of oscillators, the $++••■ 
representation of equation ( 33a | is given by 



7V0i + (N - 2) 2 



sin <pi 



CO 

K' 



(42) 



Combining (|42|) with (|33b|), the variables 0i and 2 are 

= - sin0 2 , (43) 



related through the equation 

~N<p 1 + {N-2) 2 



sm 



which admits two types of solutions 



N 



mi = 0,±1,±2, . 



2?r (2m 2 - 1) - N<f>i 
N-A 



ma = 0,1, 2, 



(44a) 
(44b) 



For each type of solution (44) there exists a function 
KZ- 1 (cf>i) which describes the fixed points of the system, 



K + l (0i, mi) = 
1 



. 47rmi \ . 47rmi 

sm 0i 1 + cos — — — ) — cos 0i sin 



N 



N 



, (45a) 



KZ 1 (4>i,m 2 ) 



sm 



(AT -2) (2m 2 - 1) 7T - AT01 
AT-4 



+ sin 0i ^ , (45b) 



where the minima d ( j >1 K± (0i) |^ 1= ^» = (maxima of 
-f^T 1 ) give the points 4>* mi and 0^ 2 where the bifurcations 
occur, 



tan ( 



cos< 



1 



cos 



An mi 
N 



N cos 



JV 

"(JV-2)(2m2-l)7r-JV</)J l 
F^4 



7V-4 



(46a) 
(46b) 



Figures |8pi and[8p show the solutions K + x (0i,mi) and 
KZ 1 (0i,m 2 ), respectively, with the stable branches in 
solid lines for N = 10. The stability of the solutions is ad- 
dressed via the eigenvalues of the 2-dimensional Jacobian 
matrix, and we see that the systems exhibits: saddle- 
(stable)node bifurcations for mi = 1 and 2, (stable)node- 
(unstable)node for mi = 0, saddle-(unstable)node for 
mi = —1 and —2 for the case KZ 1 (0i, mi). For the 
figure [8]d we obtain: saddle- (stable)node bifurcations for 
mi = and 2 on the left, saddle-(unstable)node for 
?n 2 = and 2 on the right, and a saddle-saddle m 2 = 1. 

The first unusual behavior we observe from figures [8] is 
the presence of changes in the stability of a given branch 
when the coupling varies (see mi = — 1 and —2 and the 
left branches for m 2 = and 2). From the real part 
of the eigenvalues (figure [9l we could see that A_ (0i) is 




(a) 



(b) 



FIG. 8: (Color online) Solutions of K^ 1 (4>i) representing 
the fixed points of the system synchronized region. Solid 
lines (dashed) represent stable fixed points (unstable), (a) 
Solutions of (0i,TOi) with bifurcations type node- node 
mi = 0; saddle-(stable) node bifurcation for mi = 1 and 
mi = 2, and saddle- (unstable) node for mi = —1 and 
mi = —2. (b) solutions KZ 1 (0i,m,2), with saddle-saddle bi- 
furcations m2 = 1, saddle-(stable)node bifurcation in the left 
of m,2 — and m,2 = 2 ; saddle-(unstable)node bifurcation in 
the right of m,2 = and m,2 = 2 . 



responsible for the inversions: at the bifurcations of mi = 
— 1 and —2 the fixed point is born as a saddle (Re [A + ] < 
and Re [A_] > 0), but then it changes into a stable node 
at the point where the real part of A_ becomes negative 
(figures [9^i and [9Jd) . For the bifurcations on the left of 
m 2 = and 2, shown in figure |SJd, the real part of A_ is 
positive and becomes negative after a given value of 0i, 
where the saddles become stable nodes (figures [91: and 

The exact value of the points where the branches 
change stability can be identified if we analyze all those 
solutions simultaneously as in figure |10| the branch born 
from the bifurcation on the left of m 2 = loses sta- 
bility when it touches the branch born at mi = — 1 
for K sa 0.764 with 0i = 7r/10, which becomes stable 
when the coupling increases (figure |10^ ) ; the same pro- 
cess occurs with the branch coming from the bifurcation 
on the left of m 2 = 2, which exchanges stability with the 
mi = —2 solution at K w 5.236 and 0i = — 37r/10 (figure 
Hjf)). 

Changes in the stability properties of the fixed points 
occur in the stable as well as the unstable branches as 
can be seen in the figures ^ and [9ji: at the saddle- 
unstable) node bifurcations on the right of m 2 = and 
m 2 = 2, the unstable nodes become saddles when the 
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FIG. 9: (Color online) Branches of solutions for N = 10 that 
have changed the stability by varying the strength of the cou- 
pling K. Stable solutions (unstable) are represented by solid 
lines (dashed). Upper continuous curves (purple) represents 
the real part of A_ (4>i), and down curves (blue) represent the 
real part of A+ ((/>i). Figures (a) and (b): bifurcations labeled 
by mi = — 1 and mi = — 2, respectively, showing that the 
solutions are born on the right are saddles and become stable 
node when we the real part of A_ becomes negative. Figures 
(c) and (d): bifurcations labeled by 7712 = and m,2 = 2, 
respectively, showing that the solutions are born on the left 
are stable nodes and become saddles when the real part of A_ 
becomes positive. 



(a) K~\<h) 



(b) 




ra 2 = 



ra 2 = 2 



FIG. 10: (Color online) Region where solutions exchange sta- 
bility for the cases K± x (4>i)- (a) Stable node denoted by 
m,2 = loses stability to the saddle denoted by mi = —1 
at the intersection of the solutions, (b) Stable node 7712 = 2 
loses stability to the saddle of mi = —2 in intersection of the 
solutions. Solid(dashed) lines correspond to stable(unstable) 
solutions. 



real part of A+ becomes negative. Since the solutions 
K^. 1 (tf>i) are symmetric with respect to the (f>i = n/2 
axis, the mechanism which alters the stability of those so- 
lutions is the same process as described before, with the 
exchanging points given by 4>\ — 97r/10 and <f>i = 137r/10. 

The phenomenon where stability is exchanged also 
happens on the complete phase space of the system al- 
though the exchange is not necessarily made between so- 
lutions of the same subspaces. Therefore to obtain all the 
stable solutions of the system, independent of the size of 
their basins of attraction, it is necessary to obtain all the 



points in space where two (or more) solutions collide and 
to analyze the stability of those fixed points before and 
after the collision (while at the same time the stability 
of all the other solutions should be considered). Since 
the number of solutions is exceedingly large to handle, 
even for small systems, this type of symmetry may not 
be the most enlightening choice to obtain a description 
of the synchronized region in a simple (analytical) fash- 
ion. So on the next section we will restrict our symmetry 
assumptions to more general cases, looking for the ana- 
lyticity limits of the model. 



V. SPECULAR SYMMETRY 

Let us consider an even number of oscillators with 
a configuration of frequencies that satisfy the following 
symmetry: 



^N/2+n 



n=l,..,N/2, 



(47) 



and also consider that all natural frequencies in the in- 
terval [wi,a>jv/2l are positive. If we write the variables 
4> n in terms of <Pn/2i 



sin <p n = sin ( 



K ^ ~ 3 ' 

j=n+l 



1 " 

5in<^v/2-^ 
for n = l,...,N/2 - 1 and 

1 - 

- T 



(48a) 



sm 0„ = sin i 



'N/2 



j=N/2+l 



(48b) 



for n = N/2 + 1,...,N — 1, the synchronized region is 
determined by the solutions of equation Q in the form: 



y.JV/2 



sm < 



ijV/2 ~ sm (f) N [4> N /2 ) 



K 



(49) 



From the analysis of the symmetry properties it is pos- 
sible to conclude that there is a solution of this equa- 



tion with 



i+N/2 



as long as sin < 



>N 



'AT/2- 



Since the last identity maximizes the left hand side of 



equation ( 49 ) , the critical synchronization coupling is 
given by 



N/2 



(50) 



The determination of the critical coupling with this sym- 
metry is reduced to adding the first N/2 values of the 
natural frequencies, therefore any prescription that we 
choose to determine the natural frequencies to n (N) (n — 
1, N/2) will allow us to obtain the dependence of K s 
on the size of the system N. For instance, if we con- 
sider a general case where the frequencies in the interval 
[wi, 61^/2] are obtained from a uniform distribution de- 
fined on the interval [a, b], with a > and b > a, it is 
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FIG. 11: Left: behavior of (K a (N)) for natural frequencies 
generated by a uniform distribution with a = and 6=1. 
Right: standard deviation <r [K s (N)] from the numeric out- 
come in logarithmic scale. The dots correspond to the result 
obtained through the simulation, with average obtained from 
1000 instances. The black line correspond to the curve calcu- 



FIG. 12: (Color online) Numerical solution of equation ( 52 1 in 



lated via ( 51 1 



easy to compute the mean of K s (AT) and the standard 
deviation from equation (|50|: 



(Ks(N)) 



-N, a [K s (AO] 



4VG 



N. (51) 



These results are corroborated by the simulation as illus- 
trated in figure [TT] 

Although the imposition of the symmetry defined in 



(471 on the set of natural frequencies allowed us to deter- coupling given by 



the space 4>n/2 x K representing the hxed points of the system 
with N = 6 oscillators, (a) Setting u n = ui for (wi, 0J2,^a) 
with high degeneracy, (b) The presence symmetry given by 
(471 separates the solutions such that only the bifurcations 
located in K s are degenerated. The natural frequencies are 
u)\ = 1, UJ2 = 0.7 andaj3 = 1.1. (c) The continuous symmetry 
breaking process separates the spectrum of solutions com- 
pletely, destroying the degeneracy. The natural frequencies 
are u)\ = 1, ui% = 0.7, L03 = 1.1, CJ4 = —0.8, uj 5 = — 0.75 and 
u 8 = - 1-25. 



for n = N/2 + 1, N — 1, should be considered a priori 
[22] . In order to correctly describe the systems we shall 
start with simple frequency configurations. 

The simplest configuration of frequencies satisfying the 
symmetry (47) is oj n — lo for n = 1, AT/2, with critical 



mine K s analytically, there is still the problem of finding 
all the solutions in the synchronized region. The selec- 
tion of the phase difference <Pn/2 as the independent vari- 
able warrants that no sum in the equations ( 48 1 will be 



K S (N) 



uN 



(53) 



larger than the sum on the right hand side of ( 49 1 . From 



this we infer that there exists a region of values of K 
such that the set ( 48 ) has a solution (the right hand side 
of each equation always belongs to the interval [—1,1]) 
while 4>n/2 does not satisfy (49 1. This means that all the 
solutions in the synchronized region exist after a value 
K = K s which guarantees the existence of a solution for 



If we consider a small system A^ = 6, all 2 N 2 = 16 
combinations of the kind of solutions (52) may be con- 



sidered to solve (49) numerically, as shown in figure 12 1 



equation (49 ), while in the synchronized region {K > K s ) 



the solutions of the form 4>n/2 {K) correspond to the fixed 
points of the system. 

To get all the solution in terms of <p n y^N/i) we have to 
consider the equations ( 48 ) that define relations between 



the phases, thus all possible combinations: 



Even though it is possible to determine all solutions for 
4'n/2 — 03 above K s , the high symmetry level presents 
the same problem of degeneracies in the <Pm/2 x K space 
found in the previous section: there are several values of 
phase locking (pi, 02, ^4 and 4>§ which lead to the same 
result for ^3. In larger systems the process of assigning 
an observed fixed point to a phase locked solution is prac- 
tically impossible. However there is no need for special 
prescription for the natural frequencies with in \ uj\, wjv/2] 
if we are looking for solutions that satisfy ( 47 ) . Even if 



arcsm sin < 



'N/2 



(52a) 



we consider a case with complete broken symmetry, equa- 
tion ( 49 ) will still be valid as long as the sum in the right 



hand side continues to be larger than all other sums in 
(|48>. 



7T — arcsm sin < 



'N/2 



1 " \ 

*~k £ w 4 (52b) 



j=n+l 



for n = 1, N/2 - 1, and 



Since small perturbations in the natural frequencies 
have little effect in the synchronized region, it is possible 
to start from the configuration that produces equation 
( 53 1 and continuously vary the frequencies while follow- 



( I n \ 

f>* = arcsin I sm<j) N / 2 + — ^ uij J , (52c) 

I 1 n \ 

sin0AT/ 2 + jt. w i]'( 52d ) 

\ ]=N/2+l ) 



ing the solutions in the 4>n/2 x K space. If we do this 
while preserving the specular symmetry ( 47 ) , we observe 



it — arcsm 



that all solutions born after K s are separated (as can be 
seen in figure [12)3 for u)\ = 1, 0J2 = 0.7 and W3 = 1.1). If 
we continue this process until the configurations satisfy 
only equation (491 (without symmetries), the degener- 
acy of the solutions is destroyed (the figure 12: shows 
a case with w\ = 1.0, W2 = 0.7, ^3 = 1.1, UJ4 — —0.8, 
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ljs = —0.75 and u>s = —1.25, where several solutions 
exist above K s but without any overlaps). 



( a ) 2.0 



The continuous symmetry breaking of equation ( 47 ) il- 
lustrates a property present in any configuration of natu- 
ral frequencies: the presence of degeneracies in the space 
4> n * x K is a consequence of the symmetries in the set 
{u} N . The importance of this property becomes evident 
when one consider a system with asymmetric frequency 
configuration: in general the main interest does not lie on 
the description of the full set of solutions above K s , but 
on the determination of the stable fixed points. While 
for a chain of oscillators the stable fixed point is char- 
acterized by cos0„ > (for all n), the loop structure of 
the ring allows the existence of stable fixed points with 
cos0„ < for some n (see Lee et al. [Hj), an d our 
recent work [S] showed that these phase differences may 
either be <pn* or 4>m- Hence it is only necessary to analyze 
the solutions of (49) with relations <fi n (<Pn/2) defined by 
(52a I and (52c), i.e., <j) n E [— ir,ir] for n ^ n* , while all 



other combinations of ( 52 ) always lead to unstable solu- 
tions. Now it becomes clear the choice for the $++••• in 
the previous chapter: all stable fixed points for the dy- 
namics on the whole phase space come from the solutions 
of this specific subspace. 

A general characterization of the synchronized region 
in the absence of symmetries (or with a small amount) 
was already done in our previous work jS], therefore it 
is not necessary to do it again. Nevertheless we can use 
the symmetry properties of the natural frequencies in or- 
der to study the behavior of the critical synchronization 
coupling K s as a function of the number of oscillators, as 
discussed in the next section. 



VI. ASYMPTOTIC BEHAVIOR OF THE 
CRITICAL COUPLING 



Possibly the simplest non-degenerate case analyzable 
through the solutions of ( 49 1 consists in the configuration 



where the frequencies are evenly spaced in the interval 

[-7,7]: 
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N - 1 



(-2n + JV + l), n=l,...,N. (54) 
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FIG. 13: Behavior of K s (N) obtained from different config- 
urations of natural frequencies, represented by dots, and the 
continuous line are the theoretical asymptotic curves: a) u) n is 
defined in ( 54 1 with 7=1 (logarithmic scale) and the line is 
b) Behavior of sin < 



given by equation ( 58 

c) u. 



'JN/2 - 



sin^jv ~ 2. 

is defined in (55 1 with e = 0.002 (logarithmic scale) and 
the line is given by the equation ( 62 1. d) co n is defined in ( 64 1 
with 7 = 1, and the coupling is inversely proportional to the 
system size. The curve corresponds to the asymptotic fitting 
(671 with A = 0,485. 



for n = l,...,N/2 — 1, such that the solutions in the 
synchronized region are obtained from 



sm 



» N/2 



9 AT/2 



7 iV 2 



i(N-l)K 



(57) 



The bifurcation responsible for the full synchronization 
is always the first solution of (57) and it can appear near 



any of the two boundaries of the solvability region, i.e., 
it is characterized either by sin<^)jv/2 ~ 1 or sin</>jv ~ 
— 1. The number of bifurcations giving birth to stable 
solutions is proportional to the number of oscillators, and 
for increasing N we observed that they tend to pile up 



near the minimum value of K under which equation ( 57 ) 
may present a solution. In other words, for large N it is 
possible to approximate sm0jv/2 — sin^jv sa 2 to obtain 
an asymptotic behavior for the critical coupling: 



j^-asymptotic ^jys 



7 iV 2 



8 (TV - 1)' 



(58) 



The frequencies of the oscillators are mirror images with 
respect to the axis between <f> N and 4>n/2, i- e., cj/v+i-n = Figure [T3fc shows the results obtained from the simula- 
— u) n for n = 1, N/2, which gives tion in contrast to the analytical expression in (58), while 

figure [L3|: 



'AT/2+ra 



n = l,..., N/2-1. (55) 



Assuming 4>N/2+n = 4>n the phase differences 4> n in (52) 
are given by: 



(4>n/2,K) = arcsin 



sin <j) N/ 2 - 



7 {N - 2n) 
4(N -1)K 



exhibits the actual value of the sum of the two 
sines, thus corroborating the used approximation. 

That asymptotic behavior of K s ~ N , also observed 
in the cases treated previously, is not just a consequence 
of the evenly spaced frequencies but rather from the fact 
that the interval is finite, thus the differences u„ — uj n +i 
are proportional to 1/N. If wc consider the configuration 



. (56) 



Ulr, 



(-2/1 + ^ + 1) 



1, 



,N, 



(59) 
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where u> n — uj n+ i — e, independent of TV, the equations 
(56) and (57) are written as 



»n {<PN/2, 



K) = arcsin 



sin i 



'N/2 



e(N-2nY 
8K 



for n = 1, TV/2 - 1, and 



sin 



?N/2 



2 E 0™ + Sin 



£ TV 2 



(60) 



(61) 



Since this system presents the same characteristics as the 
latter, the asymptotic limit of K s (TV) can also be ob- 



tained with sine 



'N/2 



sin (f>ff fa 2, so 



^-asymptotic ^jys 



dV 2 

^L6~ 



(62) 



as shown in figure |13b . 

For all the analyzed cases, the critical value of the cou- 
pling constant at synchronization grows with TV, which 
is common in systems with nearest neighbor interactions 
between oscillators, and this includes the chain of oscil- 
lators [IT]. In both cases given by p7| as well as the 
mirror symmetry (uiN+i-n — ~ w n) present in the evenly 
spaced configurations of frequencies ( 54 ) and ( 59 1 , any 



prescription for the distribution of frequencies that has a 
well defined behavior of the sum of the frequencies with 
TV, that is, 



N/2 

E 

n=l 



f(N), 



(63) 



will have (in general) a behavior of the critical cou pling 
of the form K s (TV) ~ \j (TV), being exact for (47) and 



asymptotic for the other cases. As in general the sum in 
(63) produces a function which increases with TV, these 
configurations will always have K s — > oo when TV — > oo. 
But, there exists at least one configuration for which this 
property is not valid: a system with mirror symmetry 
UN+i-n — —^n and natural frequencies given by 



7(-l) 



n+l 



TV- 1 



(-2n + TV + l), n = 1, N/2. (64) 



Due to the term (—1) the largest sequential sum 
J2 n w n that can be obtained corresponds to the largest 
element of the set {w} N (uii — 7), such that the equa- 
tions that describe this system in the synchronized region 
are: 



arcsin < sin 



7 [N - 2 + (-1)" (TV - 2n)} 



2(N -1)K 



for n = 2, ...,N/2, and 



20i 



N/2-1 

2 E • 

11=2 



0+0 



N/2 ' 



(65) 



1_ 
K' 

(66) 



When the number of oscillators increases the effect in 
the solvability region is the same as that in the cases of 
evenly spaced frequencies, therefore in the limit N — > 
00 we must have sin^i — sin^jv ~ 2. But as the right 



hand side of equation (66) does not depend on N, the 



critical coupling K s should decrease with TV (since for 
small systems K s > 7/2). Then, if we take the limit of a 
very large number of oscillators, we should obtain 



^asymptotic ^jy^ _ 



A 

TV 



0(TV" 



(67) 



as can be seen in figure [13p , where A is just a constant. 
The importance to obtain a configuration of frequen- 
cies which produces a behavior for the critical coupling 
K s ~ 1/TV in the LCKM resides in the comparison with 
systems that possess larger network connectivity. The 
critical coupling in the Kuramoto model (globally con- 
nected) decreases with TV for natural frequencies defined 
in a finite interval (TI5], while locally connected model 
seems to diverge in the limit TV — > 00. At first sight, 
one would expect a transition in the asymptotic limit for 
some critical network connectivity, but our result indi- 
cates that, under a proper organization of the natural 
frequencies in the LCKM, the two regimes may be con- 
nected in a continuous way, just as the type of interaction 
considered in [20]. We believe that this subject deserves 
further investigation. 



VII. CONCLUSIONS 

In the present work we gave a comprehensive descrip- 
tion of the locally coupled Kuramoto model under the 
presence of symmetries and explored to a large extent 
the analyticity limits of the system. In the following we 
summarize the salient features. 

Starting with the analysis of small systems, we gave a 
full description of the synchronized region, along with the 
determination of the critical synchronization coupling for 
some cases, and we showed how the fixed points may be 
obtained from a general configuration of natural frequen- 
cies, independent of the system size. 

To keep our attention on the analytical regime of the 
model, we considered a highly symmetric configuration, 
consisting of two natural frequencies: u> for odd num- 
bered oscillators and — w for even numbered ones. De- 
spite of its apparent simplicity, many unusual and (or) 
unexpected results were observed. 

We found that the stability of the fixed points de- 
pend on the number of dimensions of the manifold where 
the dynamics is set to happen: while for 2-dimensional 
systems (obtained from symmetric initial conditions) a 
large number of stable fixed points populate the phase 
space, the complete system is characterized by just a 
small amount of stable fixed points. 

At the onset of synchronization it is well known that 
one should expect periodic phase slips, however in the 
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two frequencies domain the system exhibits a strong de- 
pendence on the initial conditions: a regime of phase 
slips with irregular period (apparently chaotic) coexists 
with complete erratic bursting behavior characterized by 
a long time absence of phase slips as well as intermittent 
appearances. This local character of the system also ex- 
tends to the synchronized region, where we could observe 
that the stable fixed points created slightly above K s be- 
have only as local attractors in the phase space. The 
coexistence of stable fixed points may lead to the appear- 
ances of attractor crowding, as described by Wiesenfeld 
and Hadley [31], and we hope further investigation will 
clarify the issue. 

We also observed stability exchange between synchro- 
nized solutions. This feature happens when two fixed 
point collide in the phase space. The values of <f> for these 
exchanges to occur and the nature of each fixed point 
were determined through the analysis of the Jacobian's 
eigenvalues in the reduced system, and the same feature 
is expected to happen in the complete phase space. 

The presence of a type of specular symmetry enabled 
us to obtain some analytical expressions for the criti- 
cal synchronization coupling, but one important result 
came from the analysis of a continuous symmetry break- 
ing of the system: the degeneracies in the <f) n » x K space 
are lifted in the absence of symmetries on the natural 
frequencies. As a consequence we could show how just 



a small amount of the solutions will be stable, and the 
search for stable fixed points may be reduced to the anal- 
ysis of just the solutions composed by phase differences 
belonging to the interval [— tt, it]. 

By further exploring the symmetry properties we could 
manage to obtain some asymptotic results concerning 
the behavior of the critical synchronization coupling as a 
function of the system size N. While for general configu- 
ration the LCKM tends to be proportional to the number 
of oscillators, we found a counterexample where K s ac- 
tually decreases with TV. 

Many of the phenomena and results discussed here 
were not observed in a completely random natural fre- 
quencies distribution. The study of all theses symmetric 
cases have brought insights about the complexity of the 
LCKM and helped us advance in the understanding and 
the development of analytical methods to analyze this 
kind of system. 
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